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I. INTRODUCTION 

Bulk mediated surface diffusion (BMSD) defines the 
effective surface motion of particles, that intermittently 
adsorb to a surface or diffuse in the contiguous bulk 
volume. As sketched in Fig. [1] for a cylindrical sur- 
face, the particle, say, starts on the surface and diffuses 
along this surface with diffusion constant D s . Eventually 
the particle unbinds, and performs a three-dimensional 
stochastic motion in the adjacent bulk, before return- 
ing to the surface. Typically, the values of are sig- 




Figure 1: A particle diffuses in the bulk (full lines) and in- 
termittently binds to a cylinder surface on which it may also 
diffuse (broken lines). Fhis produces an effective surface mo- 
tion. We here consider the motion along the cylinder. Bottom 
right: frontal of diffusion between inner and outer cylinder. 
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nificantly larger than D s . The recurrent bulk excursions 
therefore lead to decorrelations in the effective surface 
motion of the particle, and thus to a more efficient ex- 
ploration of the surface. 

Theoretically BMSD was previously investigated for a 
planar surface in terms of scaling arguments [Jill , master 
equation schemes 3, and simulations More recently 
the first passage problem between particle unbinding and 
rebinding for a free cylindrical surface was derived [f|. 
Following our short communication [6[ we here present 
in detail an exact treatment of BMSD for a reactive cyl- 
indrical surface deriving explicit expressions for the sur- 
face occupation, the effective mean squared displacement 
(MSD) along the surface, and the returning time distri- 
bution from the bulk. In this approach different dynamic 
regimes arise naturally from the physical timescales en- 
tering our description. Thus at shorter times we derive 
the famed superdiffusive surface spreading with surface 
MSD of the form (z 2 (t)) ~ t 3 / 2 [7, 8] and the associated 
Cauchy form of the surface probability density function 
(PDF). At longer times we obtain an a priori unexpected 
leveling off of the surface MSD, representing a tradeoff 
between an increasing number of particles that escape 
into the bulk and the increasing distance on the surface 
covered in ever-longer bulk excursions for those particles 
that do return to the surface. Only when the system is 
confined by an outer cylinder eventually normal surface 
diffusion will emerge. Apart from the Levy walk-like su- 
perdiffusive regime the rich dynamic behavior found here 
are characteristic of the cylindrical geometry. 

Nuclear magnetic resonance (NMR) measurements of 
liquids in porous media are sensitive to the preferred ori- 
entation of adsorbate molecules on the local pore surface, 
such that surface diffusion on such a non-planar surface 
produces spin reorientations and remarkably long correl- 
ations times Apart from pure surface diffusion the 



2 



experiment by Stapf et al. clearly showed the influence 
of BMSD steps and the ensuing Levy walk-like superdif- 
fusion 10]. More recently NMR techniques were used to 
unravel the effective surface diffusion on cylindrical min- 
eralic rods supporting, in particular, the first passage 
behavior with its typical logarithmic dependence. BMSD 
on a cylindrical surface is also relevant for the transient 
binding of chemicals to nanotubes [ll| and for numerous 
other technological applications 0j. In a biological con- 
text, BMSD along a cylinder is intimately related to the 
diffusive dynamics underlying gene regulation (l2l . [l3| : 
DNA binding proteins diffuse not only in the bulk but in- 
termittently bind non-specifically to the DNA, approxim- 
ately a cylinder, and perform a one-dimensional motion 
along the DNA chain, as proved experimentally flil . [l5[ . 
The interplay between bulk and effective surface motion 
improves significantly the search process of the protein 
for its specific binding site on the DNA. Similarly the net 
motion of motor proteins along cytoskeletal filaments is 
also affected by bulk mediation. Namely, the motors can 
fall off the cellular tracks and then rebind to the filament 
after a bulk excursion [l6[ . Outside of biological cells the 
exchange behavior between cell surface and surrounding 
bulk is influenced by bulk excursions, the cylindrical geo- 
metry being of relevance for a large class of rod-shaped 
bacteria (bacilli) and their linear arrangements [l7j . 

The dynamics revealed by our approach may also be 
important for the quantitative understanding of coloni- 
alization processes on surfaces in aqueous environments 
when convection is negligible: suppose that bacteria 
stemming from a localized source, for instance, near to a 
submarine hot vent, start to grow on an offshore pipeline. 
From this mother colony new bacteria will be budding 
and enter the contiguous water. The Levy dust-like dis- 
tribution due to BMSD will then make sure that bacteria 
can start a new colony, that is disconnected from the 
former, and therefore give rise to a much more efficient 
spreading dynamics over the pipeline. 

In all these examples it is irrelevant which specific tra- 
jectory the particles follow in the bulk, the interesting 
part is the effective motion on the cylinder surface. We 
here analyze in detail this bulk mediated surface diffusion 
on a long cylinder. 



II. CHARACTERISTIC TIME SCALES AND 
IMPORTANT RESULTS 

In this Section we introduce the relevant time scales of 
the problem of bulk mediated surface diffusion of the cyl- 
indrical geometry presented in Fig.[T]and collect the most 
important results characteristic for the effective surface 
motion of particles. As we are interested only in the mo- 
tion along the cylinder axis z, we consider the rotation- 
ally symmetric problem with respect to the polar angle 

that will therefore not appear explicitly in the follow- 
ing expressions (compare also Sec. IIII[) . Since the full 
analytical treatment of the problem involves tedious cal- 



culations we first give an overview of the most important 
results, leaving the derivations to the forthcoming Sec- 
tions and Appendices. 



A. Characteristic time scales 

Using the result (|4"5)l for the Fourier-Laplace transform 
n(k, s) of the density of particles on the cylinder surface, 
we compute the Laplace transform of the number of sur- 
face particles, 



N a (s) 



n(z, s)dz = n(fc, s) 
No 



fc=0 



s Ai(0,s) 
D b A(0,a) 



(1) 



where k is a surface-bulk coupling constant defined be- 
low, Db is the bulk diffusion constant, 



Ai(0,s) 
A(0,s) 



Io(a$Ki(bZ) + K (aZ)h(bt), 



(2) 



and £ = ^/s/Df,. The /„ and K v denote modified Bessel 
functions. We define the Laplace and Fourier transforms 
of the surface density n(z, i) through 



and 



n(k,t) = &{n{z,t)} 



i{z,s) = 3?{n(z,t)} = 



ikz 



n(z,t)dz (3) 



t n(z,t)dt. (4) 



Here and in the following we express the transform of a 
function by explicit dependence on the Fourier or Laplace 
variable, thus, n(fc, s) is the Fourier-Laplace transform of 
n(z, t). 

From expression (JTJ we recognize that in the limit 
k — > the number of particles on the cylinder surface 
does not change, i.e., N(t) = No. The coupling para- 
meter according to Eq. (14U1) is connected to the unbinding 
time scale r fj, the bulk diffusivity Db, and the binding 
rate kb through k = Db/[kbT s\. Vanishing k therefore 
corresponds to an infinite time scale for unbinding. This 
observation allows us to introduce a characteristic coup- 
ling time 



Db 



k 2 T 2 

D h 



(5) 



Note that the binding constant kb has dimension cm/sec, 
see Section IIIII As in the governing equations the coup- 
ling constant k and the bulk diffusivity Db are the rel- 
evant parameters, the characteristic time t K in a scaling 
sense is uniquely defined. When the coupling between 
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bulk and cylinder surface is weak, k — > 0, the corres- 
ponding coupling time t K diverges. It vanishes when the 
coupling is strong, k — ¥ oo. 

While the time scale t K is characteristic of the bulk- 
surface exchange, the geometry of the problem imposes 
two additional characteristic times. Namely, the inner 
and outer cylinder radii involve the time scales 



and 



_ a 



_ b 2 



(6) 



(7) 



respectively. By definition, tb is always larger than t a . 
For times shorter than the scale t a a diffusing particle 
behaves as if it were facing a flat surface, while for times 
longer then t a it can sense the cylindrical shape of the sur- 
face. Similarly, tb defines the scale when a particle starts 
to engage with the outer cylinder and therefore senses the 
confinement. With the help of these time scales we can 
rewrite expression ([1} for the number of surface particles 
in the form 



N s (s) 



s i/2 ( s t re )i/2 + Al (0, S )/A(0 )S ) 



(8) 



where 



Ai(0, s) = Ki (V^) h (VSt&) - h (VsQ K x (Vrt~b) 

(9) 

and 
A(0,s)=J 



(Vst^j Ki (\fsh) + K (Vst^j h (\fstb) ■ 

(10) 

From the characteristic time scales t K , t a , and tb we 
can construct the three limits: 
(i) Strong coupling limit 



t K < t a < h] 



(11) 



here the shortest time scale is the coupling time. This 
regime is the most interesting as it leads to the transient 
Levy walk-like superdiffusive behavior, 
(ii) Intermediate coupling limit 



t a < t K < t b ; 



(12) 



here the superdiffusive regime is considerably shorter, 
however, an interesting transition regime is observed, 
(iii) Weak coupling limit 



t a < t b < t K . 



(13) 



To limit the scope of this paper we will not consider this 
latter case in the following. However, we note that for 
t -C tb the behavior will be similar to the t -C t K part of 
the intermediate coupling limit (ii). 



B. Important results 

We now discuss the results for the most important 
quantities characteristic of the effective surface motion. 
The dynamic quantities we consider are the number of 
particles N„(t), that are adsorbed to the inner cylinder 
surface at given time t; as well as the one-particle mean 
squared displacement 



{z\t)) 



1 

W 



z 2 n(z, t)dz. 



(14) 



This quantity is biased by the fact that an increasing 
amount of particles is leaving the surface. To balance 
for this loss and quantify the effective surface motion for 
those particles that actually move on the surface, we also 
consider the 'normalized' mean squared displacement 



\t)h 



1 



N s (t) J_ 



z 2 n(z, t)dz. 



(15) 



The detailed behavior of these quantities will be derived 
in what follows, and we will also calculate the effective 
surface concentration n(z,t) itself. Here we summarize 
the results for the surface particle number and the surface 
mean squared displacements. 



1. Strong coupling limit 

In Table U we summarize the behavior in the four rel- 
evant time regimes for the case of strong coupling. The 
evolution of the number of particles on the surface turns 
from an initially constant behavior to an inverse square 
root decay when the particles engage into surface-bulk 
exchange. At longer times, the escape of particles to the 
bulk becomes faster and follows a l/( law. Eventually 
the confinement by the outer cylinder comes into play, 
and we reach a stationary limit. 

The mean squared displacement has a very interesting 
initial anomalously diffusive behavior ~ t?l 2 0, H| . This 
superdiffusion arises due to mediation by bulk excursions 
resulting in the effective Cauchy distribution 



n(z, t) 



N Kt 



IT (z 2 + K 2 t 2 ) 



(16) 



In this initial regime we can use a simple scaling argu- 
ment to explain this superdiffusive behavior, compare the 
discussion in Ref. [l|. Thus, once detached from the sur- 
face a particle returns to the surface with a probabil- 
ity distributed according to ~ t- 1 ' 2 . Due to the diffus- 
ive coupling z 2 ~ t in the bulk the effective displace- 
ment along the cylinder is then distributed according to 
~ giving rise to a probability density ~ z~ 2 . 

Later, the mean squared displacement turns over to 
a square root behavior corresponding to subdiffusion. 
As can be seen from the associated normalized mean 
squared displacement, this behavior is due to the escap- 
ing particles. At even longer times the mean squared 
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Time regime 


N s {t) 


<* 2 (*)> 


{•2 (i))norm 


t < t K 


N 


2D s i + £> 6 * 3/2 


2D s t + — ^=D b f 3 ' /2 
3vir!i! 


t K < t < t a 


fU No 

V 7T t 1 / 2 


2D s tn + 2 ^A/ /2 


2V^D s t 1/2 + 2D b t 


t a <^t<^t b 


1 No 

2 T 


t a t K D s ^ In + VtJ^D b 


2VtJ^D s In (7^ J + 2 A,t 


t b « t 


i 6 


^D.t+^Dbt 


4"^"D s t + 2D b t 
tb 



Table I: Effective surface diffusion, strong coupling limit. For the different regimes we list the number of particles N s (t) 
on the surface, the surface mean squared displacement (z 2 (t)}, and the normalized mean squared displacement (z 2 (t)) n0 rm- 
C = exp(7) « 1.78107 where 7 is Euler's constant. 



displacement reaches a plateau value. This is a remark- 
able property of this cylindrical geometry, reflecting a 
delicate balance between decreasing particle number and 
increasing length of the bulk mediated surface translo- 
cations. This plateau is the terminal behavior when no 
outer cylinder is present. That is, even at infinite times, 
when fewer and fewer particles are on the surface, the 
surface mean squared displacement does not change. In 
presence of the outer cylinder the mean squared displace- 
ment eventually is dominated by the bulk motion and 
acquires the normal linear growth with time. 

Combining the dynamics of the number of surface 
particles and the mean squared displacement we obtain 
the behavior of the normalized mean squared displace- 
ment listed in the last column. 



2. Intermediate coupling limit 



In the intermediate coupling limit the results are lis- 
ted in Table [TTJ Also in this regime we observe the initial 
superdiffusion and associated Cauchy form of the surface 
particle concentration. The subsequent regime of inter- 
mediate times splits up into two subregimes. This subtle 
turnover will be discussed in detail below. The last two 
regimes exhibit the same behavior as the corresponding 
regimes in the strong coupling limit. 



C. Numerical evaluation 

In Figs. [5] and [3] we show results from numerical 
Laplace inversion of the exact expressions for the num- 
ber of surface particles and the surface mean squared 
displacement. We consider both the strong and inter- 
mediate coupling cases. The parameters fixing the time 
scales were chosen far apart from each other to distin- 
guish the different limiting behaviors computed in the 
following Sections. In all figures a vanishing surface dif- 
fusivity (D s = 0) is chosen for clarity. 

For strong coupling the selected time scales are t K = 
10~ 6 , t a — 1, and tf, = 10 6 in dimensionless units. There- 
fore the bulk diffusion constant becomes Df, = a 2 /t a = 25 
for our choice a = 5. The coupling constant is k = 
a/y/t a t K — 5 x 10 3 , and the outer cylinder radius be- 
comes b = a/y/tb/t a — 5 x 10 3 . 

In the intermediate regime we chose t a = 10" 6 , t K = 1, 
and tb = 10 6 . This sets the bulk diffusivity to Db — 
25 x 10 6 and the outer cylinder radius to b = 5 x 10 6 . 
These values are chosen such that we can plot the results 
for the intermediate case alongside the strong coupling 
case. 

Fig. [5] shows the time evolution of the number of sur- 
face particles, normalized to Nq = 1. For the strong 
coupling case the value remains almost constant until 
t w t K , and then turns over to an inverse square root 
decay that lasts until t « t a . Subsequently a t~ x beha- 
vior emerges. In presence of an outer cylinder, due to 
the confinement this inversely time proportional evolu- 
tion is finally terminated by a stationary plateau. In the 
intermediate coupling case similar behavior is observed, 
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Time regime 


iV s (t) 


{z 2 {t)) 


{% (^)) norm 


t < t a 


N 


2D s t+ D b t 3/2 


2L> s f+ A,i 3/2 


t a <. t < t c <IZ. t K 

t a <C t c <t < i K 


N 

transition to ~ 1/t 


2D 3 t A = 

t c ln 2 (U/[CH a ]) 

transition to ~ const 


or,,, 2A,i 2 

2L> s i H = 

t c \n 2 {U/[CH a ]) 

transition to ~ t 


t K < t < ti 


1 iV 


t a t K D s i In (j*j-J + VtaT K D b 


2Vtj:D s In \^- \ + 2£W 




— t No 

tb 


8t a t K ytatu , 

—~—D s t + 4^— D b t 


4— D s t + 2Ai 
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Table II: Effective surface diffusion, intermediate coupling limit. 



ing. The right panel of Fig. [3] shows the behavior of the 
normalized surface mean squared displacement. See Sec- 
tions [IV] and [V] for details. 



III. COUPLED DIFFUSION EQUATIONS AND 
GENERAL SOLUTION 

In this Section we state the polar symmetry of the 
problem we want to consider, and then formulate the 
starting equations for our model. The general solution 
is presented in Fourier-Laplace space. In the two sub- 
sequent Sections we calculate explicit results in various 
limiting cases, for strong and intermediate coupling. 



A. Starting equation and particle number 
conservation 

The full problem is spanned by the coordinates z meas- 
ured along the cylinder axis, the radius r measured per- 
pendicular to the z axis, and the corresponding polar 
angle 9. We are only interested in the effective displace- 
ment of particles along the cylinder axis and therefore 
eliminate the 9 dependence. This can be consistently 
done in the following way. (i) As initial condition we 
assume that initially the particles are concentrated as 
a sharp 5{z) peak on the inner cylinder surface, homo- 
geneous in the angle coordinate 9. (ii) Our boundary 
conditions are 9 independent. 

For the bulk concentration of particles in the volume 
between the inner and outer cylinders this symmetry re- 




Figure 2: Time evolution of the number of surface particles 
obtained by numerical Laplace inversion, for the following 
cases: strong binding with (green dashed-dotted line) and 
without (red full line) outer cylinder, and intermediate bind- 
ing in presence of the outer cylinder (blue dashed line). The 
indicated characteristic time scales correspond to the cases of 
strong (bottom) and intermediate (top) coupling. 



apart from the two subregimes in the range at interme- 
diate times. 

Fig.[3]depicts the behavior of the surface mean squared 
displacement. In the left panel the function (z 2 (t)} shows 
the various regimes found in the strong and intermediate 
coupling limits. Remarkably the intermediate coupling 
regime exhibits a superdiffusive behavior in the range 
t c < t < t a that is even faster than the initial i 3 / 2 seal- 
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Figure 3: Surface mean squared displacement obtained by numerical Laplace inversion. Left: The graphs show the various 
effective diffusion regimes (z 2 (t)} along the cylinder. Note the characteristic transient plateau, which is the terminal behavior 



in the absence of a confining outer cylinder. Right: Normalized function {z (t)) r _ 



The subregime t c < t < t a in the 



intermediate case distinctly shows a superdiffusive behavior that is even steeper than the initial t 3//2 scaling. We show the 
strong binding case with (red line) and without (green line) and outer cylinder, as well as the case of intermediate binding 
with outer cylinder (blue line). The characteristic time scales connected with the associated curves denote the cases of strong 
(bottom) and intermediate (top) coupling. 



quirement simply means that we can integrate out the 9 
dependence and consider this concentration as function 
of z, radius r, and time t: C — C(r,z,t). The physical 
dimension of the concentration C is [C] — 1/cm 3 . On 
the surface of the inner cylinder we measure the concen- 
tration by the density n2D(z,i), which is of dimension 
[^2d] = 1/cm 2 . Note that ri2v{z,t) does not explicitly 
depend on 9. We average this cylinder surface density 
over the polar angle, and obtain the line density n(z,t): 



n(z,t) = a / n2o(z,t)d9 — 27ran2r>(z,t), (17) 
Jo 

such that [n] = 1/cm. Note that on the inner cylinder 
with radius a the expression ad9dz is the cylindrical sur- 
face increment. The factor 2ira is important when we 
formulate the reactive boundary condition on the inner 
cylinder connecting surface line density n(z,t) and the 
volume density C(r,z,t). 

Given the line density n, the total number N s (t) of 
particles on the inner cylinder surface at given time t 
becomes 



N,(t) = 



2tv 



ad9 



ri2F>(z,t)dz = / n(z,t)dz. 



(18) 

We assume that initially No particles are concentrated in 
a 5-peak on the cylinder surface at z = 0: 



n(z, t) 



= N S(z) 



(19) 



Consequently the initial bulk concentration vanishes 
everywhere on the interval a < r < b such that 



C(r,z,t) 



= 0. 



t=o 



(20) 



Let us now specify the boundary conditions at the two 
cylinder surfaces. At the outer cylinder (r = b) we impose 
a reflecting boundary condition of the Neumann form 



^C{r,z,t) 



0. 



(21) 



In the case when we do not consider an outer cylinder 
(b — > oo) this Neumann condition may be replaced by a 
natural boundary condition of the form 



lim C(r,z,t) = 0. 



(22) 



The reactive boundary condition on the inner cylinder 
(r = a) is derived from a discrete random walk process 
in App. [A] (compare also Refs. [Hj])- Accordingly we 
balance the flux away from the inner cylinder surface, 



1 , \ 

ri2B(Z,t) 

Toff 



27rar n 



-n(z,t), 



(23) 



by the incoming flux from the bulk onto the cylinder 
surface, 



j on = lim k b C{r, z,t). 



(24) 



Here, t s with dimension [t q s] = 1/sec is the charac- 
teristic time scale for particle unbinding from the sur- 
face. It is proportional to the Arrhenius factor of the 
binding free energy e of the particles, exp(— |e|/[fesT]), 
where ksT denotes the thermal energy at temperature 
T. The binding rate kb, in contrast, has physical di- 
mension [kf,] = cm/sec, which is typical for surface-bulk 
coupling in cylindrical coordinates, compare the discus- 
sions in Refs. [H, [ID, [l8[ . For convenience, we collect the 
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coefficients in the reactive boundary condition (1241) into 
the coupling constant 



1 



2irakbT 



oil 



(25) 



such that our reactive boundary condition finally is recast 
into the form 



C(r,z,t) 



im(z, t). 



(26) 



The time evolution of the bulk density C(r, z, t) is gov- 
erned by the cylindrical diffusion equation 



C{r,z,t) = D b 



(\ d 


' d~ 


d 2 \ 


\r dr 


dr 


+ dz 2 ) 



C(r,z,t), (27) 



where D s denotes the surface diffusion coefficient. In 
many realistic cases the magnitude of D s is considerably 
smaller than the bulk diffusivity Db- The coupling term 
connects the surface density n to the bulk concentration 
C. The fact that here the bulk diffusivity occurs as coup- 
ling term stems from the continuum limit, in which the 
binding rate diverges, and therefore the binding corres- 
ponds to the step from the exchange site to the surface. 



The diffusion equations (|27j) and (|28j) together with 
the boundary conditions ([21]) and (|26|) as well as the ini- 
tial conditions (|T9|) and (|20|) completely specify our prob- 
lem. Moreover the total number of particles is conserved. 
Namely, the number of surface particles varies with time 
as 



valid on the domain a < r < b and — oo < z < oo. In 
Eq. (|27p , Di, is the bulk diffusion coefficient of dimension 
[Db] = cm 2 /sec. From a random walk perspective we 
can write Db = (<5£ 2 )/(6(<5t)), where (<5£ 2 ) is the average 
variance of individual jumps, and (St) is the typical time 
between consecutive jumps. As shown in App. [A] the 
dynamic equation for the line density n directly includes 
the incoming flux term and is given by 



d d 2 d 

—n = D s -—xn(z, t) + 27raA>— C(r, z, t) 
dt dz 1 or 



(28) 



dN a {t) 
dt 



lixaDh 



g-C(r,z,t) 



dz, 



(29) 



as can be seen from integration of Eq. (|28[) over z and 
noting that n(\z\ — > oo.t) = 0. For the number of bulk 
particles we obtain 



dN b (t) 
dt 



2nD h 



rdr 



, 1 d ( d n . 
dz -— r—C(r, z,t) 
r or V or 



f°° d 
-2iraD b J —C(r,z,t) 



dz. 



2-kDi 



r—b 



(30) 



From these two relations we see that indeed the total 
number of particles fulfills 

and therefore N s (t) + N b (t) = N . 



0, 



(31) 



B. Solution of the bulk diffusion equation 

To solve Eq. (|27l) and the corresponding boundary and 
initial value problem we use the Fourier-Laplace trans- 
form method. The dynamic equation for C(r,k,s) is the 
ordinary differential equation 

^C(r, k, s) + ~C(r, k, s) - q 2 C(r, k, s) - 0, (32) 
ar z r dr 



where we use the abbreviation 



q 2 =k 2 



D b ' 



(33) 



The reactive boundary condition becomes 



C(r,k,s) 



fj,n(k, s) 



and for the reflective condition we find 
d 



dr 



C(r, k, s) 



0. 



r—b 



(34) 



(35) 



The general solution of Eq. (|32l) is given in terms of 
the zeroth order modified Bessel functions Iq and Kq in 
the linear combination 



C(r,k,s) = AI (qr) + BK (qr). 



(36) 



The constants A and B follow from the boundary condi- 
tions, such that 



AIo(qa) + BKo(qa) = fxn(k, s) 



and 



A^-Uqr) 



+ B—K (qr) 
=b or 



= 0. 



(37) 
(38) 
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Using dIo(qr)/dr = ql\{qr) and dK Q (qr)/dr — 
—qKi(qr), we can rewrite the latter relation: 

Aqh (qb) - BqK l {qb) = 0. (39) 

The two coefficients are therefore given by 



A = 



B = ±^»n(k,s), (40) 



A(fc,s) v ' " A(fc,s) 
where we introduce the abbreviation 

A(k,s) = I (qa)K 1 (qb)+h(qb)K (qa). (41) 

Note that, due to the definition of the variable q the func- 
tion A indeed explicitly depends on the Laplace variable 
s. The solution for the bulk density C in Fourier-Laplace 
domain is therefore given by the expression 



C(r,k,s) = B^.( Kl (qb)I (qr) + A, (q h)K„ (q r) 



A(k,s) 



(42) 



C. Solution of the surface diffusion equation 

In a similar fashion we obtain the Fourier-Laplace 
transform of the dynamic equation for the surface density 
n, namely 



d 

sn(k, s) — N a — — k 2 D s n(k, s) + 2i:aD b —C(r, k, s) 

Or 

Defining the propagator of the homogeneous equation, 

1 



G s {k,s) = 



s + k 2 D s 



(43) 
>n, 

(44) 



we find 



n(k, s) = N G 8 (k, s) + G s (k, s)2itaD b —C{r, k, s) 

Or 



(45) 



From Eq. (|42|) we obtain for the reactive boundary con- 
dition that 



^-C{r,k,s) 



= _ Mn(M £^iiM (46) 



A(k,s) ' 
where 

A 1 {k,s)=K 1 (qa)I 1 (qb)-I 1 (qa)K 1 (qb). (47) 
Insertion of relation (J46J) into Eq. (j45|) produces the result 

N 



n(k, s) 



s + k 2 D s + Kq 



Ax(M 



(48) 



A(k,s) 

Here we also define the coupling constant 

K = 2wafiD b = - — — , (49) 

which allows us to distinguish the regimes of strong, in- 
termediate, and weak bulk-surface coupling used in this 
work. If we remove the outer cylinder, that enforces a fi- 
nite cross-section in the cylindrical symmetry, we obtain 
the following simplified expression, 



n(k, s) 



N Q 



k z D s 



nq 



Kxjqa) 
Ko(qa) 



(50) 



as in the limit b — > oo, we have I v (qb) — > oo and 
K v (qV) — > 0. From the Fourier-Laplace transform (|48|) 
the number of particles on the cylinder surface is given 
by 



N s (s) = n (k = 0,s), 



(51) 



following the definition of the Fourier transform. 

Plugging the result (|48|) into Eq. (j42|) we obtain the 
closed form for the Fourier-Laplace transform of the bulk 
concentration, 



C(r,fc,s) 



fiN 



A(fc, s) (s + k 2 D s + KqAi(k, s)/A(k, s)) 

Np K^Ioigr) + h{qb)Ko{qr) 

2Trak b T oS A ( fcj s ) r s + k 2 Ds + KqAl s ) / A ( fcj 8 > 



Kxiqfyloiqr) + h(qb)K {qr) 



We note that the solutions for n(k, s) and C(k, s) indeed fulfill the particle conservation, 



n{z,t)dz + 2n j dz rdr C(r, z,t) — Nq n(k,s) + 2it j rdrC(r,k,s) 

oo J a k=0 



= ^0 

k=0 S 



(52) 



(53) 



Using the results for the surface propagator n(z,t), Eq . (14"51) . we characterize the effective surface diffusion on 
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the cylinder in terms of the single-particle mean squared 
displacement 



(z\t)) =N-' / z z n{z,t)dz. 



(54) 



In Fourier-Laplace domain, we re-express this integral as 

<9 2 n(fc, s) 



-Nn 



dk 2 



(55) 



k=0 



This mean squared displacement includes the unbinding 
dynamics of particles as manifest in the quantity N s (t). 
We can exclude this effect by defining the normalized 
mean squared displacement 



<* 2 (*)>i 



N s (t) 



(At))- 



(56) 



From above results for the effective surface propag- 
ator we obtain the exact result for the surface mean 
squared displacement in App. [5] In what follows, how- 
ever, for simplicity of the argument we proceed differ- 
ently. Namely we first approximate the effective surface 
propagator n(z,t), and from the various limiting forms 
determine the surface mean squared displacement. Com- 
parison to the limits taken from the general results de- 
rived in App. [B] it can be shown that both procedures 
yield identical results. 



IV. EXPLICIT CALCULATIONS: STRONG 
COUPLING LIMIT 

In this Section we consider the strong coupling limit 
t K "C t a <C t b , representing the richest of the three 
regimes. Based on the result for the effective surface 
propagator, Eq. (14*51) . in Fourier-Laplace space obtained 
in the previous Section we now calculate the quantit- 
ies characteristic of the effective motion on the cylin- 
der surface, as mediated by transient bulk excursions. 
We consider the number of particles on the surface, the 
axial mean squared displacement, as well as the surface 
propagator. We divide the discussion into the four dif- 
ferent dynamic regimes defined by comparison of the in- 
volved time scales t K , t a , and t b . 



A. Short times, t <§C t K <§C t a <§C t b 

The short time limit t <C t K <C t a <C t b corresponds to 
the Laplace domain regime 

st K , st a , st b > 1. (57) 

1. Surface propagator in Fourier- Laplace space 

We first obtain the short time limit of the effective 
surface propagator in Fourier-Laplace space. To this end 



we note that the following inequalities hold: 



qa = ajk 2 + > aj-^- = yfsT a ^> 1, (58) 
JJb V L>b 



and thus we have 



qa 1 and qb S> 1. 



(59) 



For this case we use the following expansion of the 
Bessel functions contained in the abbreviations A(fc, s) 
and Ai(fc, s). Namely, for z — ¥ oo, 



exp(— z). 



From expressions (|4T|) and (|47p we find 

A(M~A X (M- eMq[b - a]) 



2qy ab 



(60) 



(61) 



Therefore, the surface propagator in Fourier-Laplace in 
the short time limit reduces to the simplified form 



n(k, s) 



N Q 



s + k 2 D s + Ky/k 2 + s/D b 



(62) 



2. Number of particles on the surface 

From the relation N s (s) = n{k = 0, s) we obtain the 
number of surface particles by help of the above expres- 
sion for the limiting form of n(k, s): 



N s (s) 



N Q t h 



St K + \Jst K 

Since st K 3> 1 the leading behavior follows 

s 



(63) 



(64) 



i.e., we recover that the number of particles on the sur- 
face remains approximately conserved in the short time 
regime, 



N s (t) ~ N . 



(65) 



3. Surface mean squared displacement 

The surface mean squared displacement is readily ob- 
tained from the limiting form of the surface propagator 
(l6"2l by help of relation (|55[) . Namely, we obtain 



2D S + D b / \fst~ K 



>2 ' 



(s + s/y/st^)' 

Since st K ^> I the leading behavior corresponds to 



(z 2 ( S )) ~ ^ 



(66) 



(67) 
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from which the time-dependence 



(z 2 (t)) ~ (z 2 (i)) r 



2D4 



2 D h 



3n 1 / 2 D a \t 



t 



1/2 



(68) 

yields after Laplace inversion. As in this short time re- 
gime N s (t) ~ No, the normalized surface mean squared 
displacement follows the same behavior. 

Remarkably, result (|68j) contains a contribution grow- 
ing like ~ i 3 / 2 . This superdiffusive behavior becomes 
relevant when 3> D s , which is typically observed in 
many systems. Thus for DNA binding proteins the bulk 
diffusivity may be a factor of 10 2 or more larger than the 
diffusion constant along the DNA: for Lac repressor the 
bulk diffusivity is of the order of 5..9x 10 _7 cm 2 /sec, while 
the one-dimensional diffusion constant along the DNA 
surface ranges in between 2. .9 x 10 cm /sec dHH. 



4- Surface propagator in real space 

We now turn to the functional form of the surface 
propagator n(z,t) in real space at short times in the 
strong coupling limit. We investigate this quantity in 
the limit D s = of vanishing surface diffusion. 

In the current short time limit t <C t K -C t a <C % we 
distinguish two parts of the surface density n. Let us 
start with the central part defined by k 2 s/Dj,. The 
corresponding limiting form of Eq. (|62[) is then given by 



n(k, s) 



N () 



n\k\ 



(69) 



The inverse Fourier-Laplace transform leads to the 
Cauchy probability density function 



n(z, t) 



N Kt 



tt(z 2 + kH 2 )' 



(70) 



This central part of the surfa ce p ropagator obeys the 
governing dynamic equation 0, Hflj 



(71) 



with initial condition n(z,t = 0) = Nq5(z). Here, 
we defined the space fractional derivative d/d\z\ in the 
Riesz-Weyl sense whose Fourier transform takes on the 
simple form [2~D ] 



ikz 



_d_ 

d\z 



n(z,t) I dz = —\k\n(k,t). 



Eq. (f70|) and the corresponding dynamic equation ([7T|) 
are remarkable results, which are analogous to the find- 
ings in Ref. 2] for a flat surface obtained from scaling 
arguments [22j |. It says that the bulk mediation causes 
an effective surface motion whose propagator is a Levy 
stable law of index 1. This behavior can be guessed from 



the scaling of the returning probability to the surface, to- 
gether with the diffusive scaling z 2 ~ t. However, the res- 
ulting Cauchy distribution cannot have an infinite range, 
as the particle in a finite time only diffuses a finite dis- 
tance. The question therefore arises whether there exists 
a cutoff of the Cauchy law, and of what form this is. 

The advantage of our exact treatment is that the 
Cauchy law can be derived explicitly, but especially the 
transition to other regimes studied. To this end we in- 
troduce the time-dependent length scale 



£c(t) 



(73) 



which turns out to define the range of validity of the 
Cauchy region. Namely, while at distances z > £c(t) we 
observe a cutoff of the Cauchy behavior, for z < lc(t) 
the Cauchy approximation is valid. Note that in this 
short time regime t < t K the Cauchy range scales as 
Zmax ~ y/ Dbt such that z 2 can indeed become signific- 
antly larger than n 2 t 2 at sufficiently short times, and 
thus the power law asymptotics n(z, t) ~ z~ 2 in Eq. (|70|) 
become relevant. From this Cauchy part we obtain the 
superdiffusive contribution 



z 2 nt dz 



— KyJ Dbt 
7T 



3/2 



(74) 



to the mean squared displacement, that is consistent with 
the exact forms ([68]) and (|B24|) [with D s = 0]. Cal- 
culation of the mean squared displacement however re- 
quires the k — > limit and thus involves the extreme 
wings of the distribution. As the system evolves in time 
the central Cauchy part spreads. Already in the regime 
t K < t < t a we have D^t < K 2 t 2 , and the asymptotic 



behavior 



can no longer be observed. 



To show how at very large \z\ the Cauchy form of the 
propagator is truncated we consider Eq. (|62p for small 
wave number k, 



n{k, s) 



Nps 1 / 2 

s 3/2 + X k 2 ' 



(75) 



1 /2 

where A = nD b /2. In this limit the surface propagator 
n(z,t) interestingly fulfills the time fractional diffusion 
equation 0, [H| 



33/2 



■n(z, t) 



d 2 



(76) 



with the initial conditions n(z,t = 0) = Nq5(z) and 
= 0, the second defining the initial ve- 



(72) dn(z,t)/dt 



t=o 



locity field. Here, the fractional Caputo derivative is 
defined via its Laplace transform through [25l . |26| 



QV2 

<9< 3 / ; 



■n{z, t) 



= s 3/2 n(z,s) ~ s 1/2 n(z,t = 0) 



-s-^[^-n(z,t) 



(77) 
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An equation of the form (|76| can be interpreted as a 
retarded wave (ballistic) motion [23|, [27j ■ We choose that 
the initial velocity field h(z, t)\ t _ Q vanishes. It is easy to 
show that Eq. (J7BJ) leads to the scaling (z 2 (t)) ~ i 3 / 2 of 
the surface mean squared displacement. 

The inverse Fourier transform of Eq. (|75|) leads to 



n(z, s) 



No 



exp 



,3/4 



2A l/2 s l/4 V A 1 / 2 

Inverse Laplace transform then yields 
where we use the abbreviation 



c 



AW/4 



1/4 



(78) 



(79) 



(80) 



and where M(f , /3) is the Mainardi function, defined in 
terms of its Laplace transform as [25|, [28| 



M(C,/9) = r^/ -^e'-^, 0</3<l. (81 



In the tails of the distribution, i.e., in the limit ( > 1, 
we may thus employ the asymptotic form of the Mainardi 
function, 



Mf-,/3 



a(fi) 



,(/3-l/2)/(l-/3) 



exp 



-/?) 



for 



co, where 



o(/9) 



, 6(/3) = ^/>0. 



We then arrive at the asymptotic form 



n(z,f) - Gi^^-exp 



-Co 



A 2 i 3 



(82) 
(83) 

(84) 



where Ci and C2 are positive constants. Thus, the 
Cauchy distribution in the central part is truncate d by 
compressed Gaussian tails decaying as exp (— z 4 /t 3 ) }29l |. 



B. Intermediate times t K <C t <C t a -C tb 

The range of intermediate times ( K « ( « J a « t b 
in the Laplace domain corresponds to st K <C 1 while 
st a , st b > 1. 

i. Surface propagator in Fourier-Laplace space 

As the characteristic time t K does not appear in the 
expressions A and A 1; the limiting form (|62l) is still valid 
in this regime. 



2. Number of particles on the surface 

While Eq. {63]) still holds, the leading behavior otN s (t) 
changes, as now st K -C 1: 



N s (s) ~ 



/2 : 



(85) 



and thus 



No t 1 ' 2 



In this intermediate regime the number of surface 
particles decays in a square root fashion with time t. 



3. Surface mean squared displacement 

In a similar fashion Eq. (pjijf remains valid, however, 
as we now encounter the limit st K <C 1 we obtain the 
following time dependence, 



(z 2 ( S )) 



D b tV 2 



S S 3 / 2 

After Laplace inversion the slow square root behavior 

2 



(87) 



(z*(t)) ~ 2D s t K + 



n 



1/2 



in time yields. As the number of surface particles is no 
longer constant, we obtain the normalized form of the 
surface mean squared displacement, 



2D s yM^i+2D b t 



(89) 



corrected for the square root loss of surface particles to 
the bulk, the normalized surface mean squared displace- 
ment exhibits normal diffusion. 



4- Surface propagator in real space 

In this intermediate time regime st K <C 1 and 
st a , st b 3> 1 from expression (|62p we obtain 



n(fc, s) ~ Aq- 



t 



1/2 



(90) 

Vs + D h k 2 V ; 

Recalling the translation theorem of the Laplace trans- 
form 

f{s-a)~e at f{t), (91) 
and identifying f(s) = s" 1 ' 2 , we readily find 

n(k, t) ~ Aoy^exp ( - D b k 2 t) (92) 
and thus obtain the quasi-Gaussian form 



^W^^bfc)- (93) 

This function is not normalized, corresponding to the 
time evolution of the surface particle number N s (t) ~ 
N y/t K /(TTt). 
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C. Longer times t K < t a < t < t b 

In the regime of longer times t K < t a < ^ C 4 the 
corresponding inequality in the Laplace domain reads 
st K , st a <C 1 while stb^> 1. 



1. Surface propagator in Fourier- Laplace space 

In this limit we may take qb ^> 1 and therefore have 
Ii(qb) — > 0. The ratio Ai/A is therefore approximated 
by 



Ai K^qa) 

a ~ K a { qa y 



and we find the following limiting form for the surface 
propagator, 



n(k, s) 



No 



s + k D s + nq 



K (qa) 



(95) 



2. Number of particles on the surface 
Using again the relation N s (s) = n(k = 0, s) and with 



q = ' s/ D at k = 0, we obtain 



N s (s) 



No 



(96) 



We proceed to approximate the Bessel functions in this 
expression. For small argument x, 



K {x) 



x 



ln|+ 7 



(97) 



where 7 « 0.5772 is Euler's constant. With st a ,s% -C 1 
we thus arrive at the form 



N 



(- 



(98) 



where In C = 7. After Laplace inversion (see App. [D| 
we obtain the final \ jt result for the number of surface 
particles, 



, x N , 1 



3. Surface mean squared displacement 



(99) 



The surface mean squared displacement can be ob- 
tained from expansion of the surface propagator (|95l) at 



small k. Some care has to be taken to consistently ex- 
pand the Bessel functions. We proceed as follows. Since 
qa <C 1 we make use of the expansions (j97|) and find 



n(k, s) 



No 



k 2 D s 



2 k 



a In ( 4 

111 \C>{k^a>+st a ) 



No 



k z D s 



2k 



( 94 ) n(k,s) 



We then expand in the denominator according to 

No 



k 2 D s 



2 k 



1 



a l 4 

C 2 st a 



1 



Stn 



t a In (c4l) 

(100) 

Expansion in orders of k finally leads us to 



n(fc, s) 



No 



/ Mk\ ( 4 
2 n \c? s t a 



-k 2 D J-^ln 2 . ,, , 

4 l C S*a 



2s 



.(101) 



From this expression we can now obtain the surface mean 
squared displacement in the form 



<^ 2 (s)} = - 



d 2 n(k, s) 



dk 2 



k=0 

,-. t a t K 2 ( C 2 t a 

D s In s 

2 V 4 



(102) 



Using the asymptotic Laplace transform pair (compare 
App.p 



In 2 (As) -f- - In 



Ct 



t V A 

we obtain the surface mean squared displacement 



<^(t)>~D,^ln(^- 



4t 



Dby/tatn- 



(103) 



(104) 



Normalized by the associated time evolution of the num- 
ber of surface particles the normalized surface mean 
squared displacement becomes 



4/ 



2D s Vt a t K ln ( — 



2D b t. 



(105) 



Again, this result is quite remarkable: the surface mean 
squared displacement reaches a plateau value in this re- 
gime. In absence of the outer cylinder this is the ter- 
minal behavior, reflecting the balance of ever increasing 
surface displacement due to long bulk excursions, and the 
continuing escape of surface particles to the bulk. Nor- 
malized to the time evolution of these surface particles 
we find a linear growth of the surface mean squared dis- 
placement. 
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4- Surface propagator in real space 



2. Number of particles on the surface 



In contrast to the previous two regimes, here the value 
of qa acquires values smaller and larger than 1. In the 
tails of the propagator {z ^> 1) we expand 

K (qa) ln(ga) and K x {qa) ~ l/(qa). (106) 

Therefore we can express the propagator as 

No 



n{k, s) 



(107) 



aln(l/M) 



We further approximate this expression to obtain the log- 
arithmic form 



n(k, s) ~ N -\n — = -N —\n (a 2 k 2 + st a 
k qa 2k 



(108) 



That this seemingly harsh approximation makes sense 
can be seen by evaluating the surface particle number 
N s (t) = N n(k = 0, 8 ) ~ ±NoVtat^Hl/[st a ]) leading 
to N s (t) ~ Noy/t a t K / (2t), matching our previous result 
@2). 

Formally we can now write 



n(k, t) = — ■ 



2 k 



[ e st \n(st a + k 2 a 2 



ds 
2nd 



(109) 



where the integral index Br indicates the Bromwich curve 
for the Laplace inversion. Following App. [12 we obtain 



i(k,t) = ^exp(-k 2 D b t) 



(110) 



Inverse Fourier transformation delivers the Gaussian res- 
ult 



n(z, t) 



N a a 1 



2Kt ^/ATrD b t 
with varying normalization. 



exp 



4D b t) 



(111) 



D. Long times t K <C t a <IC tb -C t 

In this final regime the outer cylinder becomes domin- 
ant, and the inequalities t K <C t a -C t b *C t correspond 
to st K ,st a ,st b <C 1 in terms of the associated Laplace 
variable. 



1. Surface propagator in Fourier-Laplace space 

In this long time regime we start with the original ex- 
pression (|48j) of the surface propagator in Fourier-Laplace 
space, and take the appropriate limits for the number of 
surface particles and the surface mean squared displace- 
ment separately. 



From Eq. (|48|) we directly obtain in the k = limit 

No 



N s {s) 



S T" K\ 



s Ai(0,g) 
D b A(0,s) 



(112) 



To calculate the approximations for Ai(0, s) and A(0, s) 
we employ the following small argument expansions of 
the Bessel functions: 



Jo (a>) » 
Ko(x) r 

K x {x)^-. 

X 



-7 — In ■ 



(113) 



Then we find 



Ai(0,s) 
A(0, s ) 



1 /i[ 

2 \ t a 
1 



lst b 



(114) 



Plugging these expansions into expression (|112[) we get 

N 



N s (s) 



and therefore 



s + s- 



tb 



t b s 



2^/tt 



N s (t) 



2N ^tJ^ 



115) 



(116) 



At long times the system reaches a stationary state due 
to the confinement by the outer cylinder. 



3. Surface mean squared displacement 

Again we start with the full surface propagator in 
Fourier-Laplace space, Eq. (l4"8"|) . and this time expand 
it around k = 0. Since st a ,st b « 1 we have qa,qb <C 1. 
To calculate Ai(fc, s) and A(k,s) at k — > we use the 
approximations (|113l) . Then, 

A l( M^#, 



A(k,s) 



1 

qb' 



(117) 



Inserting into the propagator (|48l) delivers the approxim- 
ation 

n(k,s) 1 



S + k 2 D s + ^-q 2 
Za 

1 



st b 



D s +D b 



(118) 
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We expand this expression in powers of k, obtaining 



n(k,s)~ V^_ fc 2 4^ , lK 



Stb 



2\ft a t K 



D b )+0 (k 4 ). 

(119) 

For the surface mean squared displacement we therefore 
have that 



<^ 2 (s)) = - 



d 2 n(k, s) 



dk 2 



k=0 



UatK - Ds + ADb ^k. 



and finally 



J (*)> 



17. 



(120) 



(121) 



At long times the stationary process causes a normal ef- 
fective surface diffusion. The normalized surface mean 
squared displacement attains the form 



'(*)): 



D s t + 2D b t. 



(122) 



Without surface diffusion we therefore observe normal 
linear diffusion with with the bulk diffusivity D b . In the 
presence of surface diffusion we have a correction propor- 
tional to D s . Given that t b ^> \/t a t K , the amplitude of 
this surface contribution is small. 



4- Surface propagator in real space 

In this long time regime we consider the tails of the 
propagator such that z ^> b and qa,qb <C 1. The Bessel 
functions are thus approximated by Eqs. (|113[) ; therefore, 



n{k, s) 
Since 



N 



K b 2 q 2 /(2a) 



Kb 2 



s + nb 2 (k 2 + s/D b )/(2a)' 
(123) 



2aD b 2y/UU 
we may simplify Eq. (|123[) to 



> 1 



n(fc, s) 



s + D b k 2 



(124) 



(125) 



The propagator consequently assumes the Gaussian 
shape 



n(z, t) 



1 



tb 



y/4wD b t 



exp 



4£> 6 t 



(126) 



The pref actor 2y/t a t K /t b reflects the probability that a 
certain portion of the particles is desorbed from the cyl- 
inder surface. 



V. EXPLICIT CALCULATIONS: 
INTERMEDIATE COUPLING LIMIT 

We now turn to the case of intermediate coupling 
defined by t a <C t K <C t b . 



A. Short times t < t a < t K > t b 

In this limit corresponding to st a , st K , st b < 1 we ob- 
tain the same results as in the matching limit of the 
strong coupling regime, compare Sec. IIV Al 



B. Intermediate times t a <C t <C t K -C t b 

This limit in Laplace space corresponds to the inequal- 
ities st a <C 1 and st K , st b ^$> 1. 

1. Number of particles on the surface 

In Laplace space we start from the exact expression for 
N s {s) =n(k = 0,s), 



N s (s) 



No 



"MM 

D b A(0, s) 



(127) 



where Ai(fc, s) and A(fc, s) are defined in Eqs. (j47j) and 
(|4T|) . With the asymptotic expansions of the modi- 
fied Bessel functions for small argument summarized in 
Eq. (|113p , as well as with the large argument asymptotics 
(|60| we obtain from Eq. (|127j) the result 



N s (s) 



N t c 



(128) 



ln(4/[C 2 st a }) 



where we introduced the new time scale t c = \Jt K t a which 
fulfills the inequality t a < t c < t K . 

Let us first regard the subregime t a <C t <C i c <C t K . If 
these inequalities are fulfilled, we may neglect the logar- 
ithmic term in the denominator, and find 



N s {s) 



No 

; 

S 



(129) 



i.e., the number of particles on the surface still remains 
constant to leading order: 



N s (t)~N . 



(130) 



The range t a <C t c <C t <C i K is difficult to estimate, as 
now the term linear in s and the logarithmic term in the 
denominator are of comparable order. As can be seen in 
Fig. [2j in the interval from t c to t K the number N s (t) of 
surface particles describes a quite complicated turnover 
from the persisting initial condition No to the 1 jt beha- 
vior in the following regime t K <^t <^t b . The prominent 
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shoulder visible in the double-logarithmic plot propag- 
ates to the behavior of the normalized mean squared dis- 
placement discussed below. 



2. Surface mean squared displacement 



We start from expression (|48[) for the Fourier-Laplace 
transform of the surface propagator. To determine the 
associated mean squared displacement we will need the 
small wave number approximation. Since st a <C 1 and 
stf, > lwe have qa <C 1 and qb 3> 1. With the definitions 
(|4"7)) and (|4"T]) and with the asymptotic expansions of the 
modified Bessel functions we find after a few steps 



exp(qb) f 2 



and 



(27r)Va (qb) 1/2 \ Caa 
1 exp(qb) 



Qa (2tt) 1 /2 (qb) 1/2 ' 

Thus, the following approximation 

Ai(fc,s) 1 
go In 



(131) 



(132) 



Cqa 



yields for the ratio Ai/A, and the surface propagator 
becomes 



n(k, s) 



N 



s + k 2 D s + -- 



1 



a In (2 /[Cqa]) 
At k — > we expand the logarithm as follows: 



In 



/ 2 



\Cqa 



In 



1 



Therefore 



2 VC 2 (^a + fc 2 a 2 ) 

1 / 4 \ _ fcV 

2 ^C^siJ ~ 2si a 



k 2 a 2 



In 



/ 2 



\ Cqa 



In 



C 2 si a 



(134) 



(135) 



(136) 

(133) Plugging this expansion into expression (|134|) we obtain 



n(k, s) 



This can be rephrased in the form 



2 k 



No 



a In (4/ [C 2 st a ]) 



n{k,s)=N Q 



2k 



2k a 1 



a st a In (4/[C 2 st a ]) 



2D h 



1 



Sy/tJaln (A/[C 2 st a ]) 



a ln(4/[C 2 si a ]) 



2 k 



a ln(4/[C 2 st a ]) 



(137) 



(138) 



Here we again consider the subregime t a -C t -C t c -C t K 
for which in the limit k — > 



n(fc, s) ~ iV" 



and thus 



S 



In" 



VC 2 sta 



2 2# s 4D 6 



c ^ 4 



(139) 



(140) 



After Laplace inversion we ultimately find 



(z 2 (t)) ~2D s t + 



trM' 



( 4i 



\CH a 



(141) 



According to above findings this is also the result for the 
normalized surface mean squared displacement. 

Analogous to what was said above, the following subre- 
gime t c <C t t K is difficult to estimate analytically, and 
we refer to the numerical result shown in Figs. [3J While 
for the surface mean squared displacement one can see a 
slight increase in the slope compared to the linear beha- 
vior shown by the guiding line, for the normalized analog 
we see a distinct increase in the slope after t c . 
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C. Times longer than t a <C t K 

At times longer than the characteristic scale t K our 
results again correspond to those of the strong coupling 
limit, see Sees. HVCl and HVDl 



The time evolution of P(r,z,t) is given by the diffusion 
equation 



^P(v,z,t) = D b V 2 P(r,z,t), (144) 



VI. FIRST PASSAGE STATISTICS 

In this Section we address the problem of the first pas- 
sage time statistics in our geometry, that is, the time 
it takes a particle starting at some point in between 
the two cylinders to reach the inner cylinder. As be- 
fore we neglect the dependence on the polar angle 9 
in our description. The relevant probability density is 
therefore P(r,z,t). Then 2nrP(r, z,t)drdz gives us the 
probability that the particle at time t is in the range 
(r...r + dr,z...z + dz). The initial distribution is 
smeared out on a circle of radius tq in the plane z = 0, 



P(r,z,t) 



1 



t=0 27rro 



S(r - r )5(z). 



(142) 



Here the factor l/27rro appears because of the normaliz- 
ation of the initial density, 



6 /-2ir 

rdr / d9 

Jo i X 



dz P(r, z, t) 



valid for radii a < r < b and on the entire cylinder axis, 
— oo < z < oo. The Laplace operator in polar-symmetric 
cylindrical coordinates is 



r dr \ dr I dz 2 



(145) 



In the calculation of the first passage dynamics we impose 
an absorbing boundary condition at r = a such that 



P(r,z,t) 



(146) 



while at the outer cylinder we keep the reflecting bound- 
ary condition 



d_ 

dr 



P(r,z,t) 



0. 



r—b 



(143) The result for the probability density is 



(147) 



P{r,k,s) 



1 h{qb)K Q {qr ) + K^hiqro) 
2irD b I (qa)Ki(qb) + K (qa)h(qb) 

I 



(K (qa)I (qr) - I (qa)K (qr) 



(148) 



as calculated in Appendix [C] 



A. First passage time density for times t <^.tb 

We first investigate the case when the outer cylinder 
is remote, that is, t <C t b . To this end we set b — > oo. In 
Eq. (|148[) this means that h{qb) -t oo and Ki(qb) -> 
such that 



o(qr )(^K Q (qa)I (qr) - I (qa)K Q (qr) 
2nD b K {qa) 



P(r,k,s) = 

(149) 

The probability density function for the first passage 
time is given by the radial flux through 



p(t) = 2na 



dP(r,z,t) 



dr 



dz, 



(150) 



compare also Ref. [13|. Its Laplace transform reads 

dP(r, k, s) 



p(s) = 2naD b 



dr 



(151) 



where the integral over the cylinder axis z has been re- 
placed by the zeroth Fourier mode. Inserting Eq. (I149p . 



P(s) 



qa 



+K (qa)h(qa)\ — , 

> q=^/s/D b 

K [y/st a ) I. 

+K Q (V^Q h (y/sQ } 

K {VsQ ' 



(152) 



Here we defined the diffusion time 



(153) 



r— a,fc— 



Expression (|152[) recovers a result in Ref. @. 

To evaluate this result we need the more subtle expan- 
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sion of the Bessel functions [3GJ 

K (z) ~ -7 (z){ 7 + ln(|)} + ^ 

lnz + C + 0(z 2 ) . (154) 

Here 7 ps 0.5772 is Euler's constant such that C = In 2 — 
7 > 0. We therefore find for the Laplace image of the 
first passage time density 



(s) 



ln(l/K])+2C 
ln(l/[at a ])+2C 
ln(l/[si ]) 



2C 



ln(l/[*i«]) 
ln(ioAa) 



ln(l/K]) 



1 - 



2C* 



ln(l/[st a ] 



1 



ln(l/[st a ])' 
Substituting for to we obtain 



P(s) 



1 



J Hl/[st a })' 



(155) 



(156) 



The Laplace inversion based on Tauberian theorems for 
slowly varying functions [3l| finally delivers the desired 
result 



P(*) 



, ln(r /a) 
J <ln 2 (tAa) 



(157) 



This expansion is valid in the range t > i„. We there- 
fore obtain a very subtle probability density, in which 
the logarithm ensures normalizability, however, not even 
fractional moments (t q ) with q > exist. This extremely 
shallow first passage time density is characteristic for the 
cylindrical problem. We note that in the limit tq = a we 
recover p(t) = S(t), as it should be. 



B. First passage time density for times t 3> i(> 

At times t 3> if> the outer cylinder comes into play. To 
assess the behavior of the first passage in this regime we 
insert the full solution (| 148)) into the equation (|151[) for 
the flux, finding 



P(«) 



h (yg) a (ygo) + if 1 (ygfc) Jo <ygo) 



(158) 

At st a < sto < stb < 1 we use the following expansions 
for the Bessel functions 



h 



/ s4) 

Kq (y/st a/0 ) 



1 , 

2^' 



1, 



In y/st a/0 . 



(159) 



This leads us to 



pO) 



l-st 6 ln(3t )/4 
1 - st b ln(st a )/4 



4 



1 - (t) a . 



(160) 



At long times t 3> i& wc find a finite mean first passage 
time 



(161) 



as it should be in this stationary regime. 



C. First passage time density at very short times 

t < t a 

We conclude our discussion of the first passage time 
density with the case of very short times, i < t a . In 
this regime a particle starting close to the inner cylinder 
surface does not yet feel the cylindrical geometry and we 
would naively expect that the first passage is given by 
the one-dimensional Levy-Smirnov form. 

As we may neglect the outer cylinder, we start from 
result (|149[) . We expand this result in inverse powers of 
-y/s and then perform a term-wise inverse Laplace trans- 
form. With 



we find that 
gpjVgo) 



7T , I 4^-1 / 1 

J 1 + — — + 



8z 



(162) 



1/4 



■> 1/2 
a x 

ro 

x il + 



exp 
1 / 1 



\/*0 — 

1 



st„ 



exp 



Vsto 
-ro - a 



8y/s ar 

We are thus led to the inverse Laplace transform 



(163) 



( 








a 


v. 


ro 





1/2 , 
rp - a 1 

2v^rt3 
r - a ^JTh 

exp 



exp 



ar 8v7ri 
Reorganizing this expression we find 



(ro - af 

(ro - a) 2 
4D b t 



(164) 



Pit) 



, V 2 
a \ ro — a 



ro 
x { 1 + 



V47rD 6 t 3 

#6* , 



exp 



(ro - a) 1 
4£W 



4aro 



(165) 
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Figure 4: First passage time density p(t). Left: The outer cylinder is present, eventually causing an exponential-like steep decay 
of pit), compared to the unbounded behavior. The black line shows the short time approximation (|165p . Right: Unbounded 
problem showing the long-time asymptotics. The inverse square- logarithmic term matters, as can be seen from comparison 
with the shown asymptotic behaviors. 



At very short times the first passage time density indeed 
coincides with the one-dimensional limit, reweighted by 
the ratio tJcl/tq. Note that if we keep the distance A = 
ro — a fixed but let both ro and a tend to infinity, we 
recover the result for a flat surface, 

pflatW= vAF eXP ("^)' (166) 

i.e., the well-known Levy-Smirnov distribution. 



VII. DISCUSSION 

We established an exact approach to BMSD along a re- 
active cylindrical surface revealing four distinct diffusion 
regimes. In particular our formalism provides a strin- 
gent derivation of the transient superdiffusion discussed 
earlier and explicitly quantifies the transition to other 
regimes. Notably we revealed a saturation regime for the 
MSD along the cylinder that becomes relevant at times 
above which the diffusing particle feels the curvature of 
the cylinder surface (t a ). This behavior, caused by the 
cylindrical geometry, stems from an interesting balance 
between a net flux of particles into the bulk and the 
fact that particles with a longer return time also lead 
to an increased effective surface relocation. In absence of 
an outer cylinder the saturation is terminal, while in its 
presence the MSD along the cylinder returns to a linear 
growth in time. This observation will be important in fu- 
ture models of BMSD around cylinders and particularly 
for the interpretation of experimental data obtained for 
BMSD systems. We note that in the proper limit a — > oo 
the previous results for a planar surface are recovered. 
Relaxing the strong coupling condition we demonstrated 
the existence of an almost ballistic BMSD behavior, a 



case that might be relevant for transport along thin cyl- 
inders such as DNA. 

In Ref. Q it was shown that the scaling behavior in 
the regimes below and above t a can be probed experi- 
mentally by NMR methods measuring the BMSD of wa- 
ter molecules along imogolite nanorods over three orders 
of magnitude in frequency space. For larger molecules 
such as a protein of approximate diameter 5 nm we ob- 
serve a diffusivity of 10 _6 cm 2 /sec such that for instance 
the saturation plateau around a bacillus cell (radius 1/2 
/im) sets in at around t a — 2.5 msec which might give 
rise to interesting consequences for the material exchange 
around such cells. In general, the relevance of the indi- 
vidual regimes will crucially depend on the scales of the 
surface radius and the diffusing particle (and therefore 
its diffusivity) . It was discussed previously that even the 
superdiffusive short-term behavior may become relevant 
[E 0- in general, in a given system the separation 
between the various scaling regimes may not be sharp. 
Moreover typically a single experimental technique will 
not be able to probe all regimes. It is therefore vital to 
have available a solution for the entire BMSD problem. 
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Appendix A: Derivation of the reactive boundary 
condition for a planar surface 



Similarly from Eq. (IA5|) we obtain the bulk diffusion 
equation 



We start with a derivation of the coupling between 
surface and bulk in a discrete random walk process along 
the p coordinate perpendicular to the surface, as specified 
in Fig. [5j Let Ni with i — 1,2,... denote the number of 
particles at site i of this one-dimensional lattice with spa- 
cing a. The number of particles on the surface at lattice 
site i = are termed Afo- The exchange of particles is 
possible only via nearest neighbor jumps, each charac- 
terized by the waiting time r. For the exchange between 
the surface and site i = 1 we then have the following law 



at It Tdcs 



(Al) 



where Td cs is the characteristic time for desorption from 
the surface. The bulk sites are governed by equations of 
the form 



dNx(t) 
~~dt 

dN 2 (t) 
dt 



1 

Tde; 



1 



1 



■A/o - -N x + —N 
T 2t 



2r 2t t 



(A2) 



etc. Let us define the number of "bulk" particles at the 
surface site i = through 



iVn 



2r 



■Ml- 



(A3) 



This trick will allows us to formulate the exchange equa- 
tion also for site i = 1 in a homogeneous form. Namely, 
from Eq. (|A1|) we have 



dMojt) 
dt 



1 



Moreover, from Eqs. 

dN-iit) 



2r 

2i we find 
1 



(JVi - N ) . 



and 



etc. 



dt 



dN 2 (t) 
dt 



2r 



2r 



(N - 2N t + N 2 ) . 



(Nt - 2N 2 + N 3 ) . 



(A4) 



(A5) 



(A6) 



Let us now take the continuum limit. For that pur- 
pose we make a transition from Afo — > n s as the number 
of surface particles, and Ni — > anb for the bulk concen- 
tration of particles. Expansion of the right hand side of 
Eq. (|A4I) yields the surface-bulk coupling 



dn s (t) 
dt 



a_ dn b (p,t) 
2r~ a dp 



dn b (p,t) _ J_ 2 d 2 n b 

8t ~ 2r a dp 2 ' 
P 



(A8) 



i= 2 




Figure 5: Schematic of the random walk picture of the surface- 
bulk exchange. In our derivation we pass from the discrete 
lattice with spacing a to the continuous variable p. The sur- 
face corresponds to p — 0. 



Finally, the boundary condition 



2r 



n b (p,t) 



p=0 T dc 



(A9) 



stems from our definition (IA3I). 



Appendix B: Calculation of the surface mean 
squared displacement 

To calculate the quantity we start by rewriting 
expression (1481) in the form 



n(k, s) 



N () 



s + k 2 D s + A{q, s) ' 



where 



and 



^k 2 + s/D b 



A(q,s) = nq 



Ai(g,s) 



A(g,s) 



(Bl) 



(B2) 



(B3) 



(A7) Differentiation of n(k, s) yields 



p=0 
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and thus double differentiation of n(k, s) produces 



dk 2 y ' ' ( S + k 2 D s + A(q, S )) 2 V dfc 2 



2N - 9A X L ' 



(s + fc 2 As + A(g,s)) 3 V dk 



2kD s + — . (B5) 



Since 



and therefore 
dA(q,s) 



dk 



k=0 



dq 
dk 



0, 



= and 



dn(k, s) 



Ok 



= 0, 



(B6) 



(B7) 



fc=0 



such that the first moment (z(t)} of n{z,t) vanishes, as 
it should due to symmetry reasons. Moreover, we then 
obtain 



[s + A(0, s)] 

N s( s ) 



N 2 



2D, 



d 2 A 



dk 2 



k=0- 



\ nr > d 2 A 




[ 2D ° + dk 2 


k=0- 



The second term in the square brackets can be trans- 
formed to 



d 2 A(q,s) d fdA{q,s)\dq dA(q,s)d 2 q 



dk 2 dk \ dq 
such that 

d 2 A(q,s) 



dk 2 



k=0 



dk 



Db dA(q,s) 
s 



dq dk 2 



, (B9) 



(BIO) 



k=0 



Differentiation of A(q, s) results in the expression 

q dA 1 q dA' 



dA(q,s) _ A 1 
dq K A 



1 



Ai dq A dq 



(Bll) 



(B8) At k = we have q = y/s/Db, and we ultimately obtain 
I 



N K s ) 

N 2 



2D S + k 



D b A x (0,s) 



s A(0,s) 

Here, we include the auxiliary quantities 
<9Ai 



II, - 



dq 
and 

dA 

dq 



1 3Ai 



D b VAi(0,s) dq 



1 dA 



=o A (M dq 



k=0 



(B12) 



k=0 



k=0 



K[ {y/lQ h (V^~b) - I[ (VZQ Ki (V^h) + b K x (V^Q I[ - h (V^Q K[ (Vtfb~) , (B13) 



h (V^Q Ki (V^b~) - Ki (V^Q h (V^h) + b J (y/sQ K[ (Vrt~b) + Kq (V^Q I[ (y/rtZ) , (B14) 



In these expressions the prime im plie s a derivative over 
the whole argument, and we have [30] 

I' 1 {z)=I (z)--I 1 (z), and K[(z) = -K (z) - -K^z). 

z z 

(B15) 

Taking the various limits corresponding to the cases 
discussed in Sections IIVI and |V] from above exact results 



we can indeed confirm the results obtained there based 
on the limiting forms for the surface propagator. We 
here demonstrate the corresponding derivations for the 
surface mean squared displacement based on the exact 
Eq. (|B12|) in the strong binding regime. 
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(i) At short times / <C t K <C t a <C ti, wc use the expan- 
sions 30] 



-, and K v (z) 



'2-kz ' " s ' \ 2z 

Moreover, with Eq. (|B15|> wc find that 



(B16) 



exp ( Jst a/b 

/((v^^o^)-— -jr-, ^74 (BIT) 

(27T) i /^ [St a/b ) 

(-y/sta/b) 



and 



7T 1 / 2 exp 



exp I s 



so that in this approximation 

Ai(0 )S )~A(0,s) 

thus 
dA 1 



2Va (,st Q/b ) 1/4 

(B18) 



..1/2 



f l/2 .1/2 



dq 



<9A 

9(7 



2(st a ) 1/i (st b ) 1/i 



, (B19) 



(6 -a) 



exp 



1/2 



1/2 _ 1/2 



k=0 



2(«t a )l/4( Bi6 )l/4 



and finally 



(B20) 



1 <9Ai 1 dA 



0. 



^Ai <9g A^ / k=0 
From expressions (fB12]) and ((64j) we therefore obtain 



(B21) 



/ 2/ \\ 2_D S 

With the Laplace transform pair 
1 t"- 1 



-, fc > 0, 



s k ' r(fc) 

and r(5/2) = 37r 1 / 2 /4 we arrive at the result 



(z 2 (t)) ~ 2£> s t 



1 



2JA 



3tt 1 /2 j d s U k 



1/2' 



(B22) 



(B23) 



(B24) 



where the equivalence with the normalized mean squared 
displacement is due to the fact that N s (t) ~ Ao in this 
short time regime. This result coincides with our finding 
(|68[) obtained from the approximated surface propagator. 

(ii) At intermediate times t K ; <SC t -C t a -C t& approxim- 
ation (IB2T]) still holds. From Eq. (jBT2l) and with the time 
evolution (I85[) of the number of particle on the surface in 
this regime, we obtain 



! (*)> 



N s (s) 2 
N 2 



21). n ( — 



1/2' 



i i 1/2 
2D s ^+£™, 



(B25) 



(B26) 



and, after inverse Laplace transformation, 
{ z 2 (t))~2D s t K + ^J 2 ^ 2 . 



Again, we find coincidence with the previous result 

(iii) At even longer times t K <C t a <C t <C t b we make 
use of approximations (|113j) and (|B16|) to (|B18j) . For the 
derivatives of the Bessel functions we have in this regime 



I'l (\/s*o) = Jo (Vsta) 7=7=^1 (V^T 



2' 



K[ (V^Q = -Kq (v^) - -7=^1 (y^Q l -, 



exp (js< b ] 1/2 ) 
(2^) 1 /2( s t b )i/4' 



K{ (V^h) ~ -K (V^h) 01/2A/M/4 ( B 27 ) 



2 1 /2(st fe )i/4 



Therefore, 

Ai(0,s) 

and 
9Ai 



exp ([stb] 



1/2) 



UL t k=Q 

as well as 



(2^)l/2( s t Q )l/2( s4 )l/4 

b a \ b exp ((sh) 1 / 2 ) 



A(0, S )~- 



(«t B )V2 Sia ; (2^)l/2( s4 )l/4 

1 exp([ S f 6 ]V2) 



In 



2 (2^) 1 /2( s ^)i/4 ^C 2 S t a 



(B28) 



(B29) 



(B30) 



and 

<9A 
9(7 



fc=0 



b -J 4 



exp(( s 4) 1 /2) 



(27r) 1 / 2 (si f) ) 1 /4- 
(B31) 

Here we define C = exp(7) w 1.78107, where 7 
0.577216 is Euler's 7 constant. Collecting our results 
wc find that 



1 dA l 



Ai(0,a) 9g 



1 dA 



(st a y/ 2 

and ultimately recover 



fc=0 A(O,a)0« 
-1 + - 



V 



In 



C 2 st a J J 



(B32) 



(z 2 (s)) 



N s ( s y 2 

N 2 



2D„ 



4D b 



(tj.y/ 2 ,2^4 

s m 



(B33) 
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Inserting expression (|98|) we obtain 



D s t a t K ,21 4 
— ~ — In 



C 2 st n 



D b (t a t K ) 1 / 2 



(B34) 



By Tauberian theorems the final result for the surface 
mean squared displacements in this long time regime are 



(B35) 



This result corroborates Eq. (|104l) . 

(iv) Finally, at very long times t K -C t a <C t b -C t we 
have the asymptotic behaviors 



-fo (y/st a / b ) 

K (y/st a / b ) 



1, 



-7 — In 



Va/6) 1/2, 



l '\ (\/ st a/ b ) = ^0 (y/st a / b ) - 
^1 (\J st a/b) = ~#0 (y/st a/b 



h {^/st a/b j 



1 

( S t Q/b )l/2 " ~ 2' 
^ 1 (v/s^a/fc) 



(^a/b) 1 / 2 



St a /b 



Therefore, 



0Ai 



9g 



J(^) 1/2 ln( 



C 2 st a 



and 



MM 



1 

2 V *a 



as well as 



<9A 
dq 



k=0 



— -1 



(B36) 



(B37) 



(B38) 



(B39) 



and 



A(0,s) 



For the square bracket in expression (|B12[) we obtain 



(B40) 



1 + 



2. 



and thus 



N 2 



2D S + D t 



(B41) 



(B42) 



With the help of Eq. (|115l) this implies the result for the 
surface mean squared displacement 



&t a t K 
'b 



D,t 



4(t„t«) 



1/2 



tb 



-D b t. 



(B43) 



Thus, we corroborate the finding (|121[) from the approx- 
imate calculation. With similar calculations one can re- 
produce the time evolution of the number of particles on 
the inner cylinder surface and thus the normalized sur- 
face mean squared displacement. Analogous reasoning 
confirms the results obtained for the intermediate coup- 
ling limit in Sec. [V] 



Appendix C: Calculation of the first passage time 
density 

In the Fourier-Laplace domain the diffusion equation 
(|144[) becomes an ordinary differential equation, 

^ P (r, k, s) + ~P(r, k, s) q 2 P(r, k, s) 



S{r - r ) 
2nr a D b 



where 



The boundary conditions are 
P(r,k,s) 

and 



= 0. 



(CI) 



(C2) 



(C3) 



(C4) 



r—b 



We rewrite the Fourier-Laplace transformed diffusion 
equation (|C1[) in the form 



-P(r,k,s) 



(C5) 



where the operator L and the inhomogeneity <p(r) rep- 
resent 



d 2 1 d 2 



dr 2 



and 



8{r - r ) 
2nr D b 



(C6) 



(C7) 



Eq. (|C5j) can be solved by the method of variation of 
coefficients. Namely, knowing that the solution of the 
homogeneous equation reads 



P(r,k,s)=AI (qr)+BK (qr), 



(C8) 
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to solve the full equation we assume that A and B are 
some functions of the radius r. We impose the condition 



A'{r)I (qr) + B'(r)K (qr) = 0, 



(C9) 



where the prime denotes a derivative with respect to r. 
Consequently we find 



P'(r, k, a) = A{r)I' Q {qr) + B(r)K' (qr), (CIO) 



and 



P"(r,k,s) = A{r)^{qr) + B{r)K^{qr) 

+A'(r)I^qr) + B'(r)K' (qr). (Cll) 

According to the method of variation of coefficients this 
leads to 

hP(r,k,s) =A'(r)l' Q {qr) + B'(r)K' (qr). (C12) 

With above relations we arrive at a system of two equa- 
tions with two unknowns, A' and B' , 

A'(r)I o (qr)+B'(r)K Q (qr)=0 

A'{r)I' Q {qr) + B\r)K' Q {qr) = 0(r), (C13) 

The corresponding Wronskian is [30] 

W{r) = I (qr)K' (qr) - K {qr)I' {qr) 

= -q(jo(qr)Ki(qr) + h(qr)K (qr) 

= -. (C14) 



The solutions to Eqs. (|C13|) is 



A'(r) 



W{r 



■K (qr)4>(r) 



and 



and thus 



1 



-K (qr)(j)(r)dr 



W(r)' 

' I {qr)<f>{r)dr. (C17) 



W(r) 



Mr) 
B(r) 

The general solution of Eq. (IC5[) has the form 
P(r,k,s) = 



-K (qr)<j)(r)dr + Ao I {qr) 



(C15) 



(C16) 



W(r) 



I {qr)ct>{r)dr + B K (qr), (C18) 



where the constants Aq and Bq are determined by the 
boundary conditions that are rewritten as 



P(a, k, s) = lim 



and 
dr 



P{r,k,s) 



r 1 



[ WJ?) Ko{qr ' mr ' )dr ' + Aa ) Ioiqr) + [J W(?) h{qr ' )(t>{r ' )dr ' + Ba ) Koiqr) r = 



r—b 



a W(r>) 



-K {qr'W)dr' + A ql' (qb) + 



(C19) 



^-^I^qr'Wr^dr' + B J qK' Q (qb) = 0. 
(C20) 



Since the inhomogeneity <f>{r) is a 8 distribution, in gen- 
eral we have to consider two cases separately: r < tq 
and r > tq. However, as we are interested solely in the 
first passage problem, for which we require the probabil- 
ity flow at r = a < tq we restrict ourselves to the former 
case, only. This means that Eq. (|C19|) simplifies to 



A I (qa) + B K (qa) = 0. 



(C21) 



Performing the integrals in expression (|C20[) and using 
Eq. (|C13|) we obtain the following system to determine 



the constants A$ and Bq: 



A I a (qa) + B K (qa) = 
A qh(qb) - BoqK^qb) 



q 



■(h(qb)K (qr Q ) 



2TiD b 

+K 1 {qb)h(qr )\ (C22) 



with the following determinant 



dct = -qiloiq^K^qb) + K (qa)h(qb) 



(C23) 
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The solution of Eq. (fC22|) is then 

K {qa) h(qb)K (qr ) + Ki(qb)I (qro) 



A = 



and 



Bn = 



2irD b I {qa)K 1 (qb) + K {qa)h (qb) 



Ig(qa) h(qb)K (qr ) + Ki(qb)Io(qr ) 
2nD b I (ga)lfi(g6) + Ko(qa)h(qb) 



These two relations introduced into Eq. (|C18|) we obtain 



(C24) 



(C25) 



P(r,k,s) = 



1 h(qb)Ky(qr Q ) + K^qb^grp) 
2itD b I (qa)K x (qb) + K (go) Ji (qb) 



{K (qa)h(qr)-h(qa)K Q (qr)), r<r . (C26) 



r 



This result coincides with the findings of Berg and 
Blomberg [13| , when one takes the limit of a completely 
absorbing boundary condition [in Berg and Blomberg's 
notation, this corresponds to k — ¥ oo of their reaction 
rate k]. 



Appendix D: Some Laplace transforms 

We here provide a summary of non-trivial Laplace 
transforms involving logarithmic functions, used 
throughout the text. 



1. Laplace inversion of the logarithm 

To calculate the inverse Laplace transforms in 
Eqs. d98j and ([99]). as well as (fTU5|) and (UTO]) . we use 
the direct Laplace transform of \nt. Namely, for some 
A > 0, 



A 



A 



dt 



i( 7 +ln(As)), 



(Dl) 

where 7 = — J exp(— x) \nxdx ~ 0.5772 is Euler's con- 
stant. Introducing the Laplace inversion with appropri- 
ate Bromwich path, we obtain 



Br 



qS MM ds _ set 

s 2ni \ A 



(D2) 



where C — exp(7). Differentiation of this result with 
respect to t yields 



f e st \n(As) 

JBi 



ds 
~2~Ki 



(D3) 



which is an exact result. 

Eq. (|D3]) delivers the desired result for Eq. 



S£~ x \ In 



(. 4 



1 



\C 2 st a J t' 



and for Eq. ([TTO]) : 

if- 1 {in (st a + k 2 a 2 ^j } =-i 



(D4) 



(D5) 



For the latter relation we used the shift theorem of the 
Laplace transform: 



e st \n(s + A)— = --e At . (D6) 
2ni t 



1. Laplace inversion of the squared logarithm 



To obtain the Laplace pair (|103[) we calculate 



Br 



]n(As) 



2 ds 
2-ki 



A > 0. 



(D7) 



We first deform the Bromwich path into the Hankel path 
Ha(e), that is a loop starting from —00 below the neg- 
ative real axis, merging into a small circle around the 
origin with radius |s| = e, with e — > in a positive sense, 
and finally receding to —00 moving above the negative 
real axis. The integral along the Hankel path can be 
evaluated directly: 
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'Br 
i-O 



\n(Aa] 



2 ds 
2tt7 



Ha(e) 



2 ds 
2tt7 



exp(|s|te~") \n(A\ 



s\e 



ln(As) 
2 e?|s| exp(— i7r) 



2?ri 



exp(|s|te") ln(A\s\e l7T ) 



2 rf|s| exp(z7r) 



/ e .(D8) 



Here I e — ¥ at e — )■ is the integral over the small circle To see this result let us recall the Tauberian theorems, 
of radius e around the origin. Collecting the terms we see, for instance, Ref. [Hj]. These state that if the Laplace 
find transform of some (positive) function uj(t) behaves like 



ln(j4s)— in 



2 ds 



2m 



-+c.c ■ 



so that we finally obtain 



\n(As)ds, 
(D9) 



w(s)~ s - p L - , s^O, 0<p<cx), (D13) 



Br 



ln(As) 



2 ds 



2m t 



Thus, we find Eq. (fTUgj) . 



hr 



.4 



4* 

Cta 



(D10) 

then its inverse Laplace transform has the asymptotic 
form 



7u" • (DH) 



f" 1 



(D14) 



3. Laplace inversion of the inverse square Here L(t) is a function slowly varying at infinity, i.e., 
logarithm with additional power 

We now address the Laplace pair of Eqs. (11401) and .. L(ax) /t ^-,^n 

man i mn r? r = 1 V a > 0. (D15) 

(|14ip . namely, x^oo L(x) x ' 



\C^st a J J 2 \C^t a 



[1] O. V. Bychuk and B. O'Shaughnessy, Phys. Rev. Lett. 

74, 1795 (1995). 

[2] O. V. Bychuk and B. O'Shaugnessy, J. Chem. Phys. 101, 
772 (1994). 

[3] J. A. Revelli, C. E. Budde, D. Prato, and H. S. Wio, New 

J. Phys. 7, 16 (2005). 
[4] R. Valiullin, R. Kimmich, and N. Fatkullin, Phys. Rev. 

E 56, 4371 (1997). 
[5] P. Levitz, M. Zinsmeister, P. Davidson, D. Constantin, 

and O. Poncelet, Phys. Rev. E 78, 030102(R) (2008). 
[6] A. V. Chechkin, I. M. Zaid, M. A. Lomholt, I. M. Sokolov, 

and R. Metzler, Phys. Rev. E 79 040105(R) (2009). 
[7] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000); J. 

Phys. A 37, R161 (2004). 
[8] J. Klafter and I. M. Sokolov, Phys. World 18, 29 (2005). 
[9] R. Kimmich, S. Stapf, P. Callaghan, and A. Coy, Magnet. 

Reson. Imaging 12, 339 (1994). 
[10] S. Stapf, R. Kimmich and R.-O. Seitter, Phys. Rev. Lett. 

75, 2855 (1995). 

[11] M. V. Velosoa, A.G. Souza Filhoa, J. Mendes Filhoa, and 
Solange B. Faganb, Chem. Phys. Lett. 430, 71 (2006). 

[12] P.H. von Hippel and O.G. Berg, J. Biol. Chem. 264, 675 
(1989); O. G. Berg and M. Ehrenberg, Biophys. Chem. 



15, 41 (1982); M. A. Lomholt, B. v. d. Broek, S.-M. 
J. Kalisch, G. J. L. Wuite, and R. Metzler, Proc. Natl. 
Acad. Sci. USA 106, 8204 (2009). 
[13] O. G. Berg and C. Blomberg, Biophys. Chem. 4, 367 
(1976). 

[14] I. Bonnet, A. Biebricher, P^L. Porte, C. Loverdo, O. 

Benichou, R. Voituriez, C. Escude, W. Wende, A. 

Pingoud, and P. Desbiolles, Nucl. Acids Res. 36, 4118 

(2008); I. M. Sokolov, R. Metzler, K. Pant, and M. C. 

Williams, Biophys. J. 89, 895 (2005); B. v. d. Broek, M. 

A. Lomholt, S.-M. J. Kalisch, R. Metzler, and G. J. L. 

Wuite, Proc. Natl. Acad. Sci. USA 105, 15738 (2008). 
[15] Y. M. Wang, R. H. Austin, and E. C. Cox, Phys. Rev. 

Lett. 97, 048302 (2006). 
[16] C. Bustamante, Y. R. Chemla, N. R. Forde, and D. 

Izhaky, Ann. Rev. Biochem. 73 705 (2004). 
[17] S. I. Coyne and N. H. Mendelson, Infection and Immunity 

12, 1189 (1975). 
[18] M. A. Lomholt, I. M. Zaid, and R. Metzler, Phys. Rev. 

Lett. 98, 200603 (2007); I. M. Zaid, M. A. Lomholt, and 

R. Metzler, Biophys. J. 97, 710 (2009). 
[19] R. B. Winter, O. G. Berg, and P. H. von Hippel, Biochem. 

20, 6961 (1981). 



20 



[20] 
[21] 

[22] 
[23] 
[24] 

[25] 
[26] 

[27] 



A. V. Chechkin, V. Yu. Gonchar, J. Klafter, R. Metzler, 
and L. V. Tanatarov, J. Stat. Phys. 115, 1505 (2004). 
S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional 
Integrals and Derivatives, Theory and Applications (Gor- 
don and Breach, New York, NY 1993) 
Substituting x 2 + z 2 for r 2 in Eq. (6) of Ref. @] and 
integrating x over (—00,00) produces Eq. (|70|) . 
R. Metzler and J. Klafter, Europhys. Lett. 51, 492 
(2000). 

R. Gorenflo and F. Mainardi, in A. Carpinteri and 
F. Mainardi, Editors, Fractals and Fractional Calculus 
in Continuum Mechanics (Springer- Verlag, New York, 
1997), pp. 223-276. 

I. Podlubny, Fractional differential equations (Academic 
Press, San Diego, CA, 1998). 

F. Mainardi, Y. Luchko, and G. Pagnini, Fract. Calc. 
Appl. Anal. 4, 153 (2001); Y. Luchko and R. Gorenflo, 
Fract. Calc. Appl. Anal. 1, 63 (1998); B. J. West and T. 
F. Nonnenmacher, Phys. Lett. A 278, 255 (2001). 
R. Metzler and T. F. Nonnenmacher, Chem. Phys. 284, 
67 (2002). 



[28] 



[29] 



[30] 

[31] 
[32] 



F. Mainardi, in A. Carpinteri and F. Mainardi, Editors, 
Fractals and Fractional Calculus in Continuum Mechan- 
ics (Springer- Verlag, New York, 1997), pp. 291-348. 
In our problem we deal with a truncated Cauchy distri- 
bution, therefore the normalization factor in Eq. (|70|l is 
actually different from 1/7T. However, as the truncation is 
accomplished by the rapidly decaying compressed Gaus- 
sian, the correction to 1/n is expected to be small, and 
the Cauchy part contains almost the entire probability, 



Kt , 2 t K 
—. 7r-^dz = — arctan \ — 

Z* + kH 2, TV \ t 



since t <Ci re . 

M. Abramowitz and I. Stegun, Handbook of Mathemat- 
ical Functions (Dover, New York NY, 1972) 
S. Havlin and G. H. Weiss, J. Stat. Phys. 58, 1267 (1990). 
W. Feller, An introduction to probability theory and its 
applications (John Wiley & Sons, New York, NY, 1968), 
Vol. II, Chapter XIII, Eqs. (5.15) and (5.17). 



